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For an isotropic and homogeneous nonlinear left-handed 
materials, for which the effective medium approximation is 
valid, Maxwell's equations for electric and magnetic fields lead 
naturally, within the slowly varying envelope approximation, 
to a system of coupled nonlinear Schrodinger equations. This 
system is equivalent to the well-known Manakov model that 
under certain conditions, is completely integrable and admits 
bright and dark soliton solutions. It is demonstrated that left- 
and right-handed (normal) nonlinear mediums may have com- 
pound bright and dark soliton solutions, respectively. These 
results are also supported by numerical calculations. 

Recently the study of the electromagnetic (EM) prop- 
erties of artificial complex media with simultaneously 
negative dielectric permittivity e e ff and magnetic per- 
meability l^eff has been the subject of great attention. 
Such media are usually referred to as left-handed materi- 
als (LHMs), 1 and they demonstrate a number of peculiar 
properties: reversal of Snell's law of refraction, reversal 
of the Doppler shift, counter-directed Cherenkov radia- 
tion cone, negative refraction index, the refocusing of EM 
waves from a point sourse, etc. The above-mentioned 
properties follow directly from Maxwell's equations with 
appropriate constitutive relations. Pendry 2 has proposed 
the intriguing possibility to exploit the negative refrac- 
tion index property of the LHMs in order to overcome 
known problems with common lenses to achieve a per- 
fect lens that would focus both the propagating as well 
as the evanescent spectra. 

Typical LHM are composed of a combination of a reg- 
ular array of electrically small resonant particles referred 
to as split-ring resonators (SRR's) and a regular array 
of conducting wires, 2-6 responsible for the negative fi e ff 
and e e ff, respectively. The size and spacing of the con- 
ducting elements of which the medium is composed is 
assumed to be on a scale much smaller than the wave- 
lengths in the frequency range of interest, so that the 
composite medium may be considered as a continuous 
and a homogeneous one (effective medium approxima- 
tion). Thus far, almost all properties of LHMs were 
studied in the linear regime of wave propagation, when 
both e e ff and ji e ff are considered to be independent of 
the field intensities. However, nonlinear effects in LHMs 
have been recently taken under consideration by some 
authors. 7-9 Zharov et al, 7 considered a two-dimensional 
periodic structure created by arrays of wires and SRR's 
embedded into a nonlinear dielectric, and they calculated 
e e f f and /x e / / for a Kerr-type dielectric permitivity. They 



showed that the magnetic field intensity couples to the 
magnetic resonance of the SRR in a nontrivial way, and 
that changing the material properties from left- to right- 
handed and back is allowed by varying the field intensity. 
The study of the nonlinear properties of LHMs could fa- 
cilitate future efforts in creating tunable structures where 
the field intensity changes their transmission properties. 

In the present work we show that for an isotropic, 
homogeneous, quasi-one-dimensional LHM, Maxwell's 
equations with nonlinear constitutive relations lead nat- 
urally to a system of coupled nonlinear Schrodinger 
(CNLS) equations for the envelopes of the propagating 
electric and magnetic fields. This system is equivalent to 
the Manakov model 10 that under certain conditions, ad- 
mits soliton solutions consisting of two components (vec- 
tor solitons) . For specific parameter choices, correspond- 
ing to either a left-handed or a right-handed medium, 
we find compound dark and bright soliton solutions, re- 
spectively. The constitutive relations can be generally 
written as 



D = e e// E = eE + P NL 
B = A*e//H = /xH + Mnl, 



(1) 

(2) 



where E and H are the electric and magnetic field in- 
tensities, respectively, D is the electric flux density, and 
B is the magnetic induction. The linear dielectric and 
magnetic responses of the LHM are described by e and 
H, respectively, while Pnl = cjvlE and ~M.nl — A'JViH 
are the nonlinear electric polarization and the nonlinear 
magnetization of the medium, respectively. 

It is known that e e f / and fi e f / in a LHM have to be dis- 
persive, otherwise the energy density could be negative. 1 
Their frequency dispersion, including nonlinear effects 
(but neglecting losses) , is given by 7 



e e// M = e \e D (\E\ 2 ) 



^-^onl(\H\ 2 ) 



(3) 
(4) 



where lu p is the plasma frequency, F is the filling factor, 
loqn l = u>onl(\H\ 2 ) is the nonlinear resonant SRR fre- 
quency, and £d{\E\ 2 ) = (~do + a \E\ 2 , with a the strength 
of the nonlinearity. Positive (negative) a corresponds to 
a focusing (defocusing) dielectric. For a linear dielectric 
uonl{\H\ 2 ) — > uiq, where ujq is the linear resonant SRR 
frequency. Then Eqs. (3-4) reduce to previously known 
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expressions. 5,11 The parameters F, u> p , and loq are related 
to geometrical and material parameters of the LHM com- 
ponents. Although e e ff can be readily put in the form 
e + £nl(\E\ 2 ), for /i e // this is not an obvious task, since 
fJ-eff = Veff(u 0NL ), and wqjvl depends on \H\ 2 as 7 ' 8 



£ = e(z-w't) T = e 2 t, 



(10) 



2\2 



ifl z X b \H\ z =A z E z (l-X z )(X z -fl z ) 



(5) 



where X = ujqnl/^o, Q = oj/loq, E c is a characteris- 
tic (large) electric held, and A is a function of physi- 
cal and geometrical parameters. 7 ' 8 The + (-) sign cor- 
responds to a focusing (defocusing) dielectric. Our a is 
related to the parameters of Eq. (5) as a = ±1/E 2 . 
Choosing f p — lo p /2-k = 10 GHz and f = uoo/2tt = 
1.45 GHz, left-handed behavior appears in a narrow 
frequency band (from / = 1.45 GHz to 1.87 GHz). 
For relatively small fields when fi e ff is trully held de- 
pendent, one may consider that the magnetic nonlin- 
earity jJieff = M + ^nl(\H\ 2 ) is of the Kerr type, i.e. 
Hnl{\H\ 2 ) = (3\H\ 2 (Fig. 1). The strength of the mag- 
netic nonlinearity [3 can be treated as a fitting parameter. 



=£-0.5 























u/ -0.7 






/ 


0.05 0. 1 
IHI ! /E, ! 



0.5 1 

iHr / e c 

FIG. 1. fXeff as a function of \H\ 2 /E 2 , for Q = 1.1 (solid), 
1.15 (dotted), 1.20 (dashed), 1.25 (long-dashed), A = 3, 
F = 0.4, and a > 0. Inset: Fitting to a line of the three 
first curves of the main figure for relatively small fields. 

Using Eqs. (l)-(2) and known identities, we get quite 
general vector wave equations for the fields E and H 



V 2 E - lit 



d 2 E 



(i- 



NL 



dt 2 



d .._ . TTl d 
^[(V^)xH] + - 



V(V-E) = 
Pnl) 



V 2 H - 



<9 2 H d 2 M 



— e- 



NL 



dt 2 



- V(V • H) = 



d r ,„ . d 
^[(V^)xE] + - 



e NL —(nH + M NL ) 



(6) 



(7) 



Consider an x— polarized plane wave with frequency u) 
propagating along the z— axis: 12 

E = t), 0,0), H = (0,H(z,t),0), (8) 

where 

E(z,t) = q{z,ty {k7 - ut) H{z,t)=p{z,t)e i ^ kz -° Jt \ (9) 

with k being the wavenumber. The envelopes q{z,t) and 
p(z,t) of E and H, respectively, change slowly in z and 
t. We therefore introduce the slow variables 



where e is a small parameter, and u' — doj/dk is the 
group velocity of the wave. Taking into account Eqs. (8- 
9), substituting slow variables into Eqs. (6-7), assuming 
that a — ae 2 , (5 — (3e 2 , and expressing p and q as an 
asymptotic expansion in terms of e 12 

q{€,t) = q (£,T) +sq 1 (£,T) +e 2 q 2 (£,T) H 

p(Z,r) =p (e,T)+£pi(e,r)+e 2 p 2 (e,r) + ---, (11) 

we get various equations in increasing powers of e. The 
leading order problem gives the dispersion relation u> = 
ck, where c = y/l/efi. At 0(e 1 ), the group velocity is 
given as w' = kc 2 /uj. At 0(e 2 ), we obtain 

dan ijj" d 2 qn ojc 2 , l9 a, ,9, , n 

i^ + y^ + — M?o| 2 + ef3\p \ 2 )qo = (12) 

.dp oj" d 2 p ujc 2 ~ 2 „ 2 

+ ~2 ~d£T + — W\M + «y"ko| )po = 0, (13) 

where u>" = (c? — ui' 2 )/u>, r is the slow time, and £ the 
slow space variable moving at the linear group velocity. 
By rescaling r, £ and the amplitudes q , po according to 

£ = X, T = u"t/ 2, (14) 

Q = ^\A q /aj"\q Q , P = tJ\A p /lj"\po (15) 

where A g = wc 2 /za and A p = wc 2 e/3, we get 

«QT + Qxx + (^|g| 2 + ^|P| 2 )Q = (16) 
iPT + Pxx + (a- p |P| 2 +a g |Q| 2 )P = 0, (17) 

where cr<j, p = sign(A q , p ). Eqs. (16-17) is a special case of 
the fairly general and frequently studied system of CNLS 
equations (Manakov model) known to be completely in- 
tegrable for a q = a p = a. 10 A number of bright and dark 
soliton solutions have been obtained for Eqs. (16-17) 
when a = ±1. 14-19 There is evidence that single-soliton 
(single-hump) solutions are stable while multi-hump are 
not. 16 The sign of the products fia and e/3 determine the 
type of nonlinear self-modulation (self-focusing or self- 
defocusing) effects which will occur. For a = ±1 both 
helds experience the same type of nonlinearity. 

For e, (j, > and a, (3 > we have a = +1 , and the 
system of Eqs. (16-17) accepts solutions of the form 15 



Q{X, T) = u{X)e lv o 1 P{X, T) = v{X)e lv p 



(18) 



where u, v are real functions and v q , u p are real posi- 
tive wave parameters. The latter is necessary, if we are 
interested in solitary waves that exponentially decay as 
|X| — > oo. Introducing Eq. (18) into Eqs. (16-17) we get 



uxx - v 2 u + (u 2 + v 2 )u = 
vxx - v 2 v + (v 2 + u 2 )v = 0. 



(19) 
(20) 
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FIG. 2. Envelopes qo and po of the compound bright soliton 
{<jp — cr q — +1) for v = 1 and r = A q /A v = 2 in arbitrary 
units, with the maximum amplitude of po normalized to 1. 

For v qv — u, Eqs. (19-20) have a one-parameter family 
of symmetric and single-humped soliton solutions (Fig. 

2 )15-17 

u(X) = ±v(X) = z/sech(z/X). (21) 

There are also periodic solutions of the form 

u(X) = Acos(BX) v(X) = Asin(BX), (22) 

where A = \J v 2 + B 2 with B an arbitrary parameter. 
Now without loss of generality, we take v q = \ and denote 
v p as v. This can always be acheived by a rescaling of 
variables u, v and X. Then, for < v < 1 there is 
another, in general asymmetric, one-parameter family of 
solutions for each fixed z^ 15 > 18 

u{X) = V2(l - v 2 ) cosh(i/X)/n (23) 
v(X) = -v^2{l - v 2 ) sinhpC - X )/k, (24) 

where 

k = cosh(Jf — X ) cosh(^X) — ^sinh(X — X n ) sinh(i/X) 

where Xq is an arbitrary parameter. For Xq = 0, u 
becomes symmetric and v antisymmetric. 
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FIG. 3. Envelopes qo and po of the compound dark soliton 
(<7 P = a q — —1) for k — 1 and r — A q /A p — 2 in arbitrary 
units, with the minimum amplitude of po normalized to 1. 

For e, fj, < and d,/3 > we have a = — 1 , and 
Eqs. (16-17) accept dark soliton solutions of the form 19 

Q(X, T) =P{X,T) = k [tanh(fcX) - i] e ^x-^T) _ (25) 

which are localized dips on a finite-amplitude background 
wave, as shown in Fig. 3. In this very interesting case 
of LHM the electric and magnetic fields are coupled to- 
gether forming a dark compound soliton. Note that the 



relative amplitudes are controlled by the corresponding 
nonlinearities and frequency. For a = — 1 Eqs. (16-17) 
have also solutions of the form 19 

Q(X, T) = u{X)e-< T P(X, T) = v{X)e-™l J \ (26) 

where u, v are real functions and v qi v p are real positive 
wave parameters. Introducing Eqs. (26) into Eqs. (16- 
17), we get 

u X x + v q u- (u 2 +v 2 )u = Q (27) 
vxx + v 2 p v - (v 2 + u 2 )v = 0. (28) 

For v q . v = Eqs. (27-28) allow for kink-shaped localized 
soliton solutions 20 

u(X) = ±v(X) = {v/V2) tanh(z/X/V2), (29) 

as can be seen in Fig. 4. In the context of the propagation 
of two polarization components of a transverse EM wave 
in a Kerr-type medium, the solution of this kind is often 
called a polarization domain wall. There are also periodic 
solutions of the form 

u{X) =Acos(BX) v(X) =Asm{BX), (30) 

where A = \Jv 2 — B 2 with B < v an arbitrary parame- 
ter. 
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FIG. 4. Envelopes go and po of the kink-shaped compound 
soliton (a p = a q = —1) for v = 1 and r — A q /A p — 2 in ar- 
bitrary units, with the maximum amplitude of po normalized 
to 1. 

In order to check the validity of the Kerr-type approx- 
imation for the magnetic nonlinearity, we performed ex- 
tended numerical simulations of the complete Eqs. (6 - 
7)for the fields of Eqs. (8 -9) using the full expression 
for H e ff, with ujqnl obtained analytically from Eq. (5), 
and slightly different normalization. In the simulations 
we used the exact solitons of Eq. (21) for right-handed 
medium (RHM) and Eq. (25) for LHM as initial con- 
ditions and tested their subsequent time evolution and 
resulting stability. The Q-field amplitudes are shown in 
Fig. 5 for a right- (two right figures) and a left- (two left 
figures) handed medium with £1 = 1.2 and 3.0, respec- 
tively. We take the (average) amplitude of Q to be unity 
for the right-handed medium, while the amplitude of the 
background wave of Q is unity for the LHM. Similar re- 
sults are obtained for the field P. In the left part of Fig. 
5, two bright soliton solutions of the form of Eq. (21) are 
shown, with their difference being the chosen value of v 



in the initial conditions. For relative high v the soliton 
is still stable, but its amplitude is strongly oscillating. 
For relatively low v the soliton is propagating practically 
undisturbed. In the right part of Fig. 5 two dark soli- 
ton solutions of the form of Eq. (25) are shown, with 
their difference being the chosen value of the constant 
k in the initial conditions. Again, for relatively high k 
the soliton developes strong osicllations and deforms as 
time progresses while for relatively low k the soliton is 
propagating practically undisturbed. The numerics thus 
demonstrates clearly that the analytical solutions of Eqs. 
(12-13) are good approximate solutions for the complete 
problem at relatively small fields while at larger ampli- 
tudes the exact solitons deform. This is expected since 
magnetic nonlinearity ceases to be of Kerr-type and sat- 
uration effects become more important. 




FIG. 5. Space-time evolution of solitons. Right: Bright 
solitons (RHM) with high (upper) and low (lower) amplitude 
v. Left: Dark solitons (LHM) with high (upper) and low 
(lower) amplitude k. 

In conclusion, we obtained a system of CNLS equa- 
tions equivalent to the Manakov model describing the 
propagation of EM waves in a nonlinear LHM (as well as 
in a nonlinear RHM) for relatively small fields. Unlike 
a recent article, 21 where only magnetic nonlinearity was 
considered in the propagation of EM waves in LHMs, in 
the present work both magnetic and dielectric nonlincar- 
ities were retained. The present analysis however does 
not address the issue of the switching effect between nor- 
mal medium (RHM) and LHM properties. Although not 
obvious from the Eqs. (3 - 5), the reduction to the Man- 
akov form becomes possible after a proper approximation 
of the complex Eqs. (4) and (5). It turns out that for the 



choice of parameters corresponding to LHM, the system 
admits compound dark soliton solutions. For the choice 
of parameters corresponding to normal medium, the sys- 
tem admits compound bright soliton solutions. Refer- 
ence to compound " dark" and " bright" solitons we mean 
that both soliton components, i.e. both the envelopes of 
the electric and the magnetic fields, are either both dark 
or both bright. The described effects are not limited 
to the specific SRR - wire system. They should also be 
present in other nonlinear LHM designs, such as photonic 
crystals, 22 coupled nanowire systems, 23 the transmission 
line systems, 24 and photonic systems. 25 The case where 
the fields in the medium experience different type of non- 
linearity, leading to <r q = —o v = 1 or a v = —a q = 1 cor- 
responds to a medium of positive e e / / and negative [i e f / 
or negative e e / / and positive /i e / / , respectively. This in- 
teresting case where the Manakov system does not seem 
to be integrable will be treated numerically in a following 
publication. 
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